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Abstract 

Electromagnetic modes with parabolic-cylindrical symmetry and their dynamical variables are 
studied both in the classical and quantum realm. As a result, a new dynamical constant for the 
electromagnetic field is identified and linked to the symmetry operator which supports it. 
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I. INTRODUCTION 



The dynamical variables of the electromagnetic (EM) field define its mechanical iden- 
tity and are essential for understanding the effects of the field on charged particles. As a 
consequence, finding the natural dynamical variables of the EM field and their relations to 
mode properties has been directly linked to the development of classical and quantum EM 
theory. Historically, the photon concept emerged from suggesting a definite relationship 
between the energy (linear momentum) of a photon and the frequency uj (wave vector k) of 
plane waves. Similarly, the relationship between the angular momentum of EM waves and 
their polarization is very important for understanding, e. g., atomic processes mediated by 
photons. 

Photons associated to EM modes with non Cartesian symmetries are characterized by 
sets of dynamical constants different from those of plane waves. An example corresponds to 
circular-cylindrical EM waves known as Bessel modes l|] or their paraxial analogue, i. e., 
Laguerre- Gaussian beams j^. Bessel photons carry a well defined orbital angular momen- 
tum 2, 3] proportional to the winding number m of their vortices {4, 5|. Another example 
corresponds to Mathieu modes which exhibit elliptical-cylindrical symmetry. Mathieu pho- 
tons carry constant values for the balanced composition of the orbital angular momentum 
with respect to the foci of the elliptical coordinate system jol . 

The purpose of this paper is to analyze the mechanical properties of the EM modes for the 
last coordinate system with translational symmetry along an axis known to have separable 
analytical solutions, i. e., the parabolic-cylindrical coordinate system. The solutions of 
the corresponding wave equation can be expressed in terms of Weber functions, giving this 
name to the EM modes. We show that a balanced composition of a component of the lineal 
momentum with a component of the angular momentum is a natural dynamical variable for 
these modes. Weber photons carry a well defined value of this variable. 

Weber beams of zero-order have already been experimentally generated by means of a 
thin annular slit modulated by the proper angular spectra [3]. This setup was conceived as 
a variation of that originally used by Durnin et al. l|] for generating Bessel beams. Higher 
order Weber beams can also be produced by holograms encoded on plates 0] or in spatial 
light modulators 8|. 

Parabolic-cylindrical coordinates. 
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The parabolic-cylindrical coordinate system {u, v, z) is defined by the transformations 9| 

X + i y = -{u + i vY , z = z (1) 

where x, y and z are the well known Cartesian coordinates, and u G (—00, 00) and v G [0, 00). 
Surfaces of constant u form half confocal parabolic cylinders that open towards the negative 
X axis, while the surfaces of constant v form confocal parabolic cylinders that open in the 
opposite direction. The foci of all these parabolic cylinders are located at x=0 and y=0 for 



each z value. The scaling factors associated to u and v are hu = = h = \/v? + v"^. In 
the following, the notation represents the unitary vector related to a given coordinate x, 
the shorthand notation is used for partial derivatives with respect to the variable x, and 
do =■ dct = \dt with c the velocity of light in vacuum. 

II. A PARABOLIC SCALAR FIELD AND ITS DYNAMICAL VARIABLES 

The scalar wave equation has separable solutions invariant under axial propagation, 

V^^ = dl^, *(f, t) = 7/;(rl)e'('=^"-"*Vv^ (2) 

in four coordinate systems: Cartesian, circular-, elliptic- and parabolic-cylindrical coordi- 
nates. For parabolic-cylindrical symmetry Helmholtz equation reads 

[h-'{dl + d',)+kl]ij{u,v) = 0, kl = e-kl (3) 

where the real constants k = u/cis the magnitude of the wave vector for a given frequency 
uj, kz its axial component and k_[_ its perpendicular component. If il){u,v) = U{u)V{v) 

{dl + klu^ - 2k±a) U{u) = 0, (4) 
{dl + klv^ + 2k^a) V{v) = 0, (5) 

with 2kj_a the separation constant. These equations are known as parabolic cylinder or 
Weber differential equations. Solutions for this differential set can be expressed as Frobenius 
series, parabolic cylinder functions, Whittaker functions, Hermite functions, and others 
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12] . Here, the solutions are expressed in terms of confluent hypergeometric functions 



of the first kind,iFi: 



V,,uUv) = s/f^e-^"^ + ^; ^Q, (7) 

with (u = k±u'^, = k±v'^, Ue = 1 {uo = 3) for even (odd) parity functions: Ue{—u) = Ue{u) 
{Uo{—u) = —Uo{u)). The normahzation factors are taken as 

-y/ TT sec(za7r) ■\/27rsec(ia7r) 

So 



^ |r(3/4-za/2)|' " |r(l/4-za/2)r 
These expressions reduce directly to the Frobenius series, presented in Ref. [10| and used 



in Ref. 



13l | to introduce parabohc optical wave fields in the paraxial regime. In order to 



guarantee that these scalar fields vanish when the absolute values of the coordinate variables 



tend to infinity, a must be real 



The set {^p,K^ = {kz, uj, a)} is complete and orthogonal. 



Each function "ifp^^ satisfies the eigenvalue equations, 

= i-ir^p,^{u,v,z,t) (8) 

-idz '^p^^{u,v,z,t) = kz'^p^^{u,v,z,t) (9) 

idt '^p,^{u,v,z,t) = uj'<iJp^^{u,v,z,t) (10) 
2A ^p^,{u,v,z,t) = ^—dl-—d^^p^^{u,v,z,t) 

= 2k±a <I/p,^{u,v,z,t) (11) 

The operator A is directly identified as a generator of the balanced composition of a rotation 
around the z-axis and translations along the y-axis since 

A = (1/2)9,. + ydly - xdl = [hpy + pyh) II (12) 

where Iz = —i{r x V)z, and py = —idy. For scalar fields and spacetime continuous sym- 
metries, the generators of infinitesimal transformations turn out to be good realizations of 
the corresponding dynamical operator. In that sense, A can be related to the product of 
the z-component of the angular momentum and the y-component of the linear momentum. 
Then, the eigenvalue equation Eq. ffTTl) means that the scalar field '^p^kz,uj,a carries a well 
defined value of that momenta product. 

Travelling scalar Weber modes are defined by 



^Mr,t) = (^e,.(w,^) ±#o,K(w,^))(e^^'=^^""*Vv^, (13) 
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a) b] 




X [units of A] X [units of A] 



Figure 1: (Color Online) Sample of (a) transverse intensity and (b) polarization structure of an 
odd EM TE Weber field. They correspond to the eigenvalue a = —2, kz = 0.995a;/c and the unit 
of length is taken as the wave length. 

and these modes are orthonormal |3| 

/ n,Am±A^ = '^Kh - K)6{k± - k'JS{a - a'). 



III. PARABOLIC-CYLINDRICAL EM MODES 

In Coulomb gauge, any given solution for the vector electromagnetic field, A, can be 
written as a superposition of modes related to a complete set of solutions for the scalar 
wave equation, \[^, which are identified as Hertz potentials For EM fields in parabolic- 
cylindrical coordinates having a well defined behavior under 

Ip,, = 4^^)M^p,, + (14) 

with the vector operators given by the expressions 

The constants A^if^ and A^k^^ are proportional to the amplitude of the transverse electric 
(TE) and transverse magnetic (TM) EM fields, as can be directly seen from their connection 
with the associated electric and magnetic fields, E = —ddA and i? = V x A yielding 

4,. = Af^^d^m,,, - 4™)9,,Mvi/,,,. (16) 
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Similar expressions can be written for the travelling EM modes associated to Eq. (fT3i) . The 
intensity and polarization structure of a EM Weber beam are illustrated in Fig. [1] The 
instantaneous electric field orientation has a nontrivial structure and, as a function of time, 
E preserves its direction while its magnitude oscillates at each point. 



IV. DYNAMICAL CONSTANTS FOR PARABOLIC-CYLINDRICAL EM MODES 

Given a symmetry generator Noether theorem is usually applied to the field Lagrangian 
density 

1 ^ 

= — V (9^ A, - d^A^) {d^A^ - d^A^) , (17) 

in order to find a dynamical constant for the electromagnetic field. This theorem states that 
jisl l , if under an infinitesimal transformation of the space coordinates — > + ^'^^^p 
and the field ^4^ — > + (^P^dojp, the Lagrangian is left invariant, then the current 

K (18) 

AO- 

has a null divergence. As a consequence, 0° defines the density of a dynamical variable 
whose integrated value over a volume can change only due to the flux of the current 
through the surface that delimits the volume. 

Since we are working with EM modes that have a particular symmetry the subset of 
all possible transformations whose generators are directly identified from the eigenvalue 
equations (IHl fTTl) is particularly relevant. Under an infinitesimal translation along the main 
direction of propagation Sz or a time translation 6t, the EM field changes according to the 
expressions A^ ^ A^ + d^A^Sz or A^ ^ A^ + dtA^St, respectively. This is reflected in the 
fact that the field momentum-like variable 

Pi^'^'-'^' = ^ 1 d'x{E^l X ^ = TE, TM, (21) 

is independent of time if the integration is taken over the whole space. Similarly, the energy- 
like integral 

= ^jd'^ [4:1 ■ 4'!.' + 4:^ ■ 4t'] (22) 



is also constant. In fact, pj^'^'^^'P '"^ ^ and £!(*'P>'''P'''^') are proportional to each other with k^/u) 
the constant of proportionality. Notice that, in both cases, the dynamical constant can be 
inferred from the factor defined in Eq. (fT9|) up to a term proportional to the divergence 
of a vector field, for instance 

^? = E ^^^^^^ = ^ - E ^^(^^^^•)- (23) 

i i 

contains just the transformation of the field A^. 

Under an infinitesimal rotation around the z-axis with an angle 600, the electromagnetic 
field A has a well defined transformation rule, Ai = Ai ^ eij^AjSu, which is independent of 
the origin of space coordinates. If one considers the Noether term Eq. (fT9l) . an expression 
for the helicity is found: 

Si^'^''^'^'^''^'^ = (I/Attc) [ d'x (E« X A^;:i,),. (24) 

Jv 

This is another dynamical constant for Weber EM modes as can be directly verified by 
substituting the general expression for the vectors E and A in terms of the Hertz parabolic 
modes. As expected, it turns out that for a given mode p, k, the helicity Sz is different 
from zero only if the amplitudes and A^"^^^ that define the polarization of a mode 

are complex. The Noether term A°, Eq.( !20l) . defines the density of the 2;-component of the 
orbital angular momentum 



For Weber modes, similarly to periodic plane waves 16|, = V ■ G + 6Cz with the latter 
term becoming zero for p = p', k, = n' . Explicitly, 

^= (26) 

i 

^ j=u,v 

Thus, a not null orbital angular momentum for a given Weber mode {p, k) in a volume V 
can be due just to a flux of the vector field G through the boundary surface. 
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Let us consider A as generator of a transformation for the vector field A. Noether theorem 
as described above concerns first order differential operators as generators of continuous 
symmetries while A contains second order operators. Nevertheless, the S° term associated 
to this transformation gives rise to 

E(4:^)^-^(^S«')^- = ■ ^p'^^' + ^ ■ ^^^^ (27) 

3 

The first resulting term in Eq. (1271) is proportional to the integrand that defines the energy, 
Eq. fl22|) . Thus, for integrations over a finite volume V, C*-*-* defines the flux of 21 through 
the surface around the integration volume. Eq. (l27j) supports the identification of 21 as the 
electromagnetic dynamical variable related to the generator A. 

As for the discrete symmetry, using the properties of the scalar function U under the 
reflection of u and the expression of the EM modes in terms of Hertz potentials, it is 
straightforward to find the reflection properties of the electric field Ep^^ for each mode: 




V. QUANTIZATION OF THE EM FIELD IN TERMS OF WEBER MODES 

Standard quantization rules require a proper normalization of the EM modes so that each 
travelling photon carries an energy hw. The classical electric field amplitude is substituted 
by the electric field per photon times the creation operator for the given travelling mode, 
Ak^ = ef^ak\ I^Kp = h/k'^k'j_. The quantum energy and the momentum along z operators 
take the form: 

S = J2^Ni^\ '^-^ = Y. ^^^^ (29) 
in terms of the number operator: 

A^f = \ {o^'^ + 0^^") , W^.o^'] = kA,^', (30) 
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allowing the identification of hkz and hio with the photon momentum along z and the photon 
energy respectively. As for the helicity, 

= E ^ (^(Ti.)t^(TM) _ )ar)t) . (31) 

K 

A quantum analysis of the relation between polarization and helicity can be carried out in 
analogy with the study in Ref. for Bessel fields. Finally, in the quantum realm the field 
operator associated to A is 

^ = J2h^k^aNi'\ (32) 

An overall factor h was introduced so that the dynamical variable 2t for a photon has units of 
linear momentum times angular momentum as expected for the quantum variable associated 
to pyh- 



VI. DISCUSSION 



The parabolic-cylindrical modes differ from other separable cylindrical modes by hav- 
ing {IzPy +Pylz)/'2 as a symmetry operator. We showed that the quantum numbers of the 
EM modes {kz,uj,a} are related to their linear momentum along z, the energy, and the 
symmetrized product of the angular momentum along z and the momentum along y. The 
helicity, a property intrinsic to the vector nature of the EM field, was shown to be diago- 
nal in the circular basis resulting from the complex superposition of TE and TM modes. 
The dynamical variable 2t is gauge dependent although it can be written in a gauge in- 
dependent looking form for monochromatic modes for which Ep^^ = iujAp^^. Contrary to 
standard EM dynamical variables which depend on products of E and i?, 21 depends on 
the products of E, B, and their derivatives. Since Weber EM modes form a complete set, 
21 = (l/47rc) d^xEjAAj will be a conserved quantity for any EM wave A whenever the 

flux of 21 through asymptotic parabolic cylinder surfaces at infinity is null. This flux can 
be evaluated writing the given EM field A as a superposition of Weber modes and applying 
Eq. ([27D. 

The mechanical effects of Weber beams on cold atoms deserve a detailed study both 
quantum mechanically and semiclassically (in complete analogy to that already done for 
Mathieu beams j^). However, there are some qualitative features that can be expected 
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without performing such an analysis. For instance, since under stationary conditions, cold 
noninteracting atoms in a red-detuned light beam have a higher probability of being located 
in the higher intensity regions of the beam, the corresponding squared atomic wave function 
mimics the intensity pattern of the light field. Thus, for a Weber lattice, the atomic wave 
function is expected to have a geometrical structure similar to that of the scalar Weber 
function, Eq. ([2]). This structure gives rise to the eigenvalue equations Eqs. (IHl fTT]) . Thus, 
necessarily 21 define a natural dynamical variable for the mechanical description of the atomic 
cloud in a Weber lattice. A careful analysis concerning this idea is in progress. 
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